arXiv:astro-ph/9706117vl 12 Jun 1997 


A&A manuscript no. 


(will be inserted by hand later) 


ASTRONOMY 

AND 

ASTROPHYSICS 

1.2.2008 

Your thesaurus codes are: 

07(08.01.3; 08.13.1; 08.13.2; 08.18.1; 08.05.1) 



Unstable waves in winds of magnetic massive stars 

Henning Seemann and Peter L. Biermann 

Max-Planck-Institut fiir Radioastronomie, Auf dem Hiigel 69, 53121 Bonn, Germany 
Received November 25, 1996; accepted June 2, 1997 


Abstract. We use a luminous fast magnetic rotator model to 
analyze the influence of a magnetic field on the linear waves 
induced in the wind of a massive star by the radiative instability. 
We show that a twisted magnetic field can drive a strong wind 
with a wind efficiency {Mvoo)/{L/c) > 1 even without mul¬ 
tiple scattering. The radiation amplified waves in the wind are 
modified by the twisted magnetic field so that they can enhance 
the wind and lead to overestimates for M and Voo- Finally we 
argue that the spin down time might be consistent with the life¬ 
time derived from the mass loss rate within the uncertainties 
regarding the stellar structure. Therefore our model may help 
to explain high, observed values for M and Voo without being 
ruled out by the spin down problem. 

Key words: Stars: atmospheres - Stars: magnetic fields - Stars: 
mass-loss - Stars: rotation - Stars: early-type 


1. Introduction 


The wind of massive stars is primarily driven by the inte ractio n 
with stellar radiation in many lines (Lucy & Solomon 1970 ). 
Based on this concept, the theory of Castor, Abbott, and Klein 
(1975 hereafter CAK) can explain many aspects of hot star 
winds. But the basic model of CAK fails to describe the strong 
winds observed in some massive stars, including Wolf-Rayet 


stars. Beside multiple scattering effects (Gayley et al. 1995) a 


rotating magnetic field as driving force were proposed (Friend 


& MacGregor 1984 Poe et al. 1989, Maheswaran & Cassinelli 


1992, and Biermann & Cassinelli 1993) to overcome this prob¬ 


lem. 


We emphasize, that the observations of nonthermal radio 


emission from OB (Bieging et al. 19891 and Wolf-Rayet stars 


(Abbott et al. 1986) give clear evidence for non neglible mag¬ 


netic fields on these stars. On the other hand direct observations 
give only rather high limits on the magnetic field (Landstreet 


1982), especially for very early-type stars. Landstreet’s sample 


contains only one O-star (09.5), which is not massive enough to 
form a Wolf-Rayet star later. Further indirect hints for a strong 
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magnetic field come from radio o bserva tions of superno va ex- 
plosions ( Bierm ann & Cassinelli 1993 Biermann et al. 1995i 
Biermann 1997). More observations (e.g. Ignace et al. 1995) 


are necessary, before our model can be verified or ruled out. 

Previous models for rotating magnetic stars suffer from 
one serious problem: A sufficient enhancement of the wind 
is connected to a severe loss of angular momentum. The star 
would spin down faster than is a llowe d from stellar evolution. 
But Maheswaran & Cassinelli ( 1994 ) argued that a massive 
star could even reach the Wolf-Rayet phase as a fast magnetic 
rotator. Magnetic fields may also play an im portan t role in 
rotationally compressed winds (Cassinelli et al. 1995). 

Another important feature of hot star winds is the instability 
of the radiation force, which leads to a highly perturbed wind. 
The stability and propagation of linear hydrodynamical waves 
in the ra diatio n field of the star were analyzed by Owocki & 
Rybicki ( |l984| herea fter OR , and 1985 ), who u nified prior work 
by MacGregor et al. ( |1979| ) and Abbott ( 1980 ). They found that 
the linear waves grow very rapidly and therefore cannot be neg¬ 
lected in an overall model. The resulting waves can influence 
the stellar wind in two ways: (i) They can gain momentum and 
energy from the radiation field, transport it, and later dissipate 
it into t he win d. This can help to drive the unperturbed wind 
(Koninx 1992) e.g. by increasing the real M at the base of the 
wind, where shocks might not have formed yet. (ii) Waves can 
also steepen into shocks and influence the observation of fun¬ 
damental wind parameters: The observed terminal velocity Voo 
is inferred from the blue edge of P-Cygni lines. This maximal 
velocity has to be compared with the unperturbed Uoo plus the 
velocity amplitude of the waves or, if outward running waves 
have already steepened into outward running shocks, to the un¬ 
perturbed Voo plus the shock speed. The last scenario is more 
probable due to the high amplification rates found in OR. Krolik 
& Raymond (1985) and MacFarlane & Cassinelli ( 1989| ) ana¬ 
lyzed the structure of a single outward running shock shell in a 
nonmagnetic wind and found that the shock shell is driven by 
the radiation field to a velocity much higher than the velocity of 
sound. The observed values for the mass loss rate M are influ¬ 
enced by waves in two ways: (i) M is inferred from observed 
values of the density phy M ^ pr^Voo- An overestimated Voc 
therefore leads to an overestimated M. (ii) The values for p are 
inferred from observed radio or UV fluxes. A matter distribution 
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perturbed by waves or shocks generates higher fluxes, which 
lead to higher values for p if the perturbations are not taken into 


account properly (Abbott et al. 1981 and Hillier 19911. Shocks 


are also a good source for the observed X-rays from hot stars. 


Lucy (1982) developed a heuristic model consisting of a chain 


of outward running forward shocks mainly in order to de scribe 
the X-ray emission from hot stars. But Owocki et al. ( 1988 i 
found in a nonlinear time-dependent calculation, that waves in 
the stellar wind are dominantly running inward and steepen 
into inward running, reverse shocks. In this case the arguments 
about Voc and also Koninx’s model do not apply, so that such 
waves can not significantly help to explain the high values for 
M and Voo observed in the wind of massive stars. 

The aim of this paper is to show in the l imit o f a linear anal¬ 
ysis that the objections of Owocki et al. (1988) do not apply, 
if a strong magnetic field and rotation are present. Therefore 
we join the ideas of a rotating magnetic held (Maheswaran & 
Cassinelli 1992,1994) and of unstable waves (Owocki & Ry- 


cous, perfectly conducting fluid as described in Jackson (1962 1 


cally thick winds. We start with 
dp 

0 = ^ + V • {pv) (1) 

dv B 

0 = p-^ + p{v ■ V)t) -I- Vp -I- — X (V X B) - F,ad (2) 


dB 

0 = —— V X {v X B) 


(3) 


where Brad is the force term due to radiation pressure. Now we 
express the pressure p by c^p, where Cj is the speed of sound, 
and replace p, v, and B by sums of a equilibrium value and a 
perturbation. In the comoving reference frame {v = 0)13 this is; 


bicki 1984 )OR)) in the wind of a massive star and do a linear 


stability analysis. We show that the magnetic held increases 
the number of wave modes and changes their properties signif¬ 
icantly. High phase velocities can be achieved due to the high 
Alfven velocity. Therefore inward running waves will not be 
advected outward and will not steepen into reverse shocks any¬ 
more, which dominate at large radii. Furthermore outward and 
inward running waves have the same growth timescale in the 
short wavelength regime where both modes are unstable. So we 
can expect outward running waves far from the star and forward 
shocks in the nonlinear regime. This model does not e xclude 
multiple scatteri ng as described e.g. by Lucy & Abbott ( 1993 i 
or Gayley et al. ( 1995 ). Rather both could be joined to create a 
model for an even stronger wind. 

In Sect. 2 we derive the dispersion relation for radiatively 
amplified waves in the presence of a magnetic held. In Sect. 3 
we discuss our unperturbed wind models. In Sect. 4 we discuss 
the waves we found for the wind models of Sect. 3. In Sect. 5 
we discuss the observational consequences of our model. And 
in Sect. 6 we describe some conclusions. 


2. The wave equations 

To analyze waves in the wind of hot stars we start from the equa¬ 
tions of magnetohydrodynamics for a compressible, nonvis- 


p = p + 5p[x,t) 

V — Sv{x,t) 

B = B + SB{x,t) 

Equations 1-3 give to first order in small quantities: 
dSp{x, t) 


(4) 

(5) 

( 6 ) 


0 = 


dt 


+ pV ■ 5v{x, t) 


0 = + cl'S/5p{x, t) + ^ X (V X 5B{x, t)) 


0 = 


dSB{x, t) 

Wt 


— V X {Sv{x,t) X B) 


(7) 

(8) 
(9) 


We assume now, that the radiative force acts only in radial di¬ 
rection and depends only on the radial velocity. These equations 
can be reduced to a single equation for Sv: 


0 = 


d^5v{x,t) 2 


SB 


— Cj V(V • 5v{x, t)) + 


Da X V X [V X {Sv{x,t) X Da)] 
1 ddVr{x, t) SFi-ad 
p dt dVr 


( 10 ) 


and add the analytic description of OR for the influence of 
the stellar radiation on the plasma. OR was criticized by Lucy 


(1984) for the neglect of damping due to diffuse radiation. But 
Owocki & Rybicki ( 1985| ) showed that this effect reduces the 
amplification rate only by approximately 50% at one stellar ra¬ 
dius and 20% at infinity. In this initial paper we emphasize the 
effect of the magnetic field. Therefore we choose the simple de¬ 
scription of OR instead of the more exact b ut also more involved 
description of Owo cki & Rybicki ( 1985 ) or the description of 
Gayley & Owocki (1995), who analyzed the instability in opti- 


with the vectorial Alfven velocity Da = B/^/Anp. If we use the 
result of OR for the linear perturbation of the radiative force, 
we can get a dispersion relation from this wave equation for 
plane waves: 

Sv{x,t)=6ve''^-^-'^‘^*. (11) 

Equation 10 then becomes: 

0 = —cu^Sv + (Cj -I- v^)(k ■ Sv)k + (Da ■ fc)[(DA • k)dv 

Slkr 

— {va • 6v)k — (k ■ i5u)ua] — ^— —uj{6v ■ er)er (12) 

X “b 

Here we used the analytic form of the perturbation of the radi¬ 
ation force derived in OR. 


n = 


1 — 1/2 


1 — aCAK f^Sobo 


( 13 ) 


* Variables marked with a refer to the unperturbed wind. 
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is their combined amplification rate of all lines. 


3. The unperturbed wind models 


X = 


/2c“cak-1/2 I 

1 — acAK .^Sobo 


a 

Vr 


(14) 


is their mean blue-edge absorption strength with the empirical 
parameter c « 1.6. And is the radial velocity of the unper¬ 
turbed wind in the rest frame of the star. Lsobo = t^th/ {dvr/dr) 
is the Sobolev length. Equation 12 is a vector equation, which 
is linear in 6v. We can think of it as a generalized eigenvalue 
problem: 

(A(/c, ua, Cj ) — a;B(fcr, ii, x) — = 0 (15) 


We can find lo and 5v numerically. Then SB and 5p follow 
from Eqs. 7 & 9: 

5p = —{k ■ 6v) (16) 

iJ 

SB = -Uk ■ Sv)B - (B ■ k)Sv] (17) 

iJ 

Although Eq. H is very involved in the general case, we can 
find an analytical solution for a simplified situation: 


k{r) = ker (18) 

VA{r) = vab^ (19) 


The latter approximation is quite accurate far away from the 
star. In this limit we find 


Q,k 


± 


Q,k 


■i\) +(Cs+t^A)^^- (20) 

2(x + ifc) VV2(x + ifc)y' 

In the long wavelength limit {k x), where the waves are 
stable, this leads to 


CO = 


2 


k. 


( 21 ) 


In the case of a weak magnetic field with Ur S> Cs ^ va, this 


resembles Abbott’s (1980) result for stable radiative-acoustic 


waves with a fast inward and a slow outward mode. In the case 
of a strong magnetic field we have Vr ^ va ^ Cs- This leads to 
higher phase velocities for both modes and reduces the relative 
difference between inward and outward waves. Additionally 
inward running waves are not advected outward by the average 
wind motion any more, because their phase velocity is higher 
than the velocity of the unperturbed wind. In the short wave 
length limit (k ^ x) we find 


.£i , 

CO = 1— ± 
2 



+ (C? -b 


( 22 ) 


These waves propagate, if Q. is less than 2k^Jcl -b v\, with the 
same phase velocity inward and outward. The amplification rate 
is also the same for both modes. 

In the case of no magnetic field Eq. [pleads to 


0 = w^-b 


Q.kr 


ikr 


— dk^co — 


Q,k. 


X + ik, 


+ kl). (23) 


Eor kg = k^ = 0 this reproduces the result for isothermal waves 
found in OR. Eor oblique waves there is a third wave mode. 


To analyze the effect of these waves we construct three wind 
models for a standard massive star. Table [| gives the fundamen¬ 
tal parameters of our standard star. The first model (model A) 


Table 1. Fundamental parameters of our standard star 


M 

23Mq 

L 

1.7 X lO^L© 

R 

8.5i?0 

T 

60000K 

OCAK 

0.56 

kcAK 

0.28 


is the standard analytic CAK wind for a star without mag¬ 
netic field, rotation, or pressure. In this model we reproduce 
the previous result of do minantly inward running waves found 
by Owocki et al. ( 1988 ). In the second model (model B) we 
add a radial magnetic field. Since this magnetic field is paral¬ 
lel to the natural stream lines of the wind of a nonrotating star, 
this field does not change the velocity profile of the unperturbed 
wind. But it changes the microphysics for waves. The last model 
(model C) is a luminous fast magnetic rotator model. The rotat¬ 
ing magnetic field provides an additional driving force, which 
changes the properties of the wind drastically. The terminal 
velocity decreases and the mass loss rate increases. The wind 
efficiency {Mvoo) / {L / c) is 3.8, whichis much higher than for 
a purely radiatively driven wind in the single scattering limit, 
where the efficiency can not be higher than unity. In spite of the 
known spin down problem we choose a strong magnetic field 
and a high rotation rate in order to emphasize the influence of 
these parameters. Eor this model we use the equations derived 
by Biermann & Cassinelli (1993), which do not assume that 
the magnetic field is radial close to the star. This leads to a 
smaller and more accurate value for the Alfvenic radius Ra, 
which controls the angular momentum loss of the star. Table ^ 
gives the magnetic field, rotation rate, and the resulting values 
for the above mentioned unperturbed wind models. In Sect. 5 
we will discuss how these values and their observation can be 
influenced by waves. In this paper we do not discuss the in¬ 
fluence of the waves back on the unperturbed wind as Koninx 
( 1992 ) did. A model with rotation but without magnetic held is 
not included, because this model has the same microphysics for 
waves as model A. Just M and the velocity dependence on r are 
different. Model C differs from model B in the local conditions 
by the fact that, due to the rotational twist, the magnetic field 
and the direction of the radiative force are not parallel anymore. 
Even at the base of the wind is approximately — 1 .61 Br- 
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Table 2. Results for the unperturbed wind models 


Model 

A 

B 

C 

BrO [G] 

0 

500 

500 

ri/ i^crit 

0 

0 

0.92 

Woo [km s“'] 

1146 

1146 

784 

M [10“'’M0 yr-‘] 

0.6 

0.6 

17 

-Rt— 2/3 [-^o] 

8.5 

8.5 

11 

{Mv^)/{L/c) 

0.2 

0.2 

3.8 


4. Numerical results for waves 

Since our unperturbed wind model is limited to the equatorial 
plane we limit our discussion for waves to the same plane. 
We have still two free parameters then: The wavelength and 
the azimuthal angle of k. In the first section we discuss radial 
waves. This will show all important properties of this wave 
model. In the second section we discuss briefly the influence of 
the azimuth angle of k. 


4.1. Radial Waves 


Figures |1|&Q show nimerical results for wind model A&B. 
The sonic wave modesEl found for model A are identical to the 
slow magnetosonic modes of model B. For model B we And 
additionally four fast wave modes. These modes, the fast mag¬ 
netosonic and the Alfvenic modes, have the same phase velocity 
uph ~ Da, because both, B and fe, are parallel to e^. But these 
wave modes are stable and show no dependence on wavelength. 
Therefore we concentrate our discussion for model A&B on the 
slow waves: Fig. shows the dependence of the waves on the 
wavelength for r = 2i?, plotted relative to the Sobolev length, 
which is the relevant length scale for radiative wave amplifi¬ 
cation in the model of OR, we use here. R — 8.5i?0 is the 
stellar radius. In the long wavelength limit we And stable waves 
with a high phase velocity inward and a low phase velocity 
outward. This reproduces Abbott’s ( 1980 ) result of radiative- 
acoustic waves. In the short wavelength limit we And the same 
amplification timescales and phase velocities (except the di¬ 
rection) for both modes. In this limit we would expect to see 
both wave modes in the wind. But the most interesting case 
is the bridging case A « Tsobo- Here we And the shortest am¬ 
plification timescale for all wavelengths. Inward waves are two 
orders of magnitude faster in amplification than outward waves. 
They have also a higher phase velocity. This resembles the re¬ 
sult of Owocki et al. (1988), who did a nonlinear calculation 
and found that primarily the inward running waves steepen into 
reverse shocks, which are advected outward. The aim of this 


^ Ordered by the radial phase velocity we denote in our figures the 
fast magnetosonic modes with 1&6, the Alfvenic modes with 2&5, and 
the slow magnetosonic modes with 3&4. The sonic modes of model A, 
which are identical to the slow magnetosonic modes of model B, are 
denoted with 3&4 as well. 


paper is to argue that this scenario changes if a magnetic field 
and rotation are involved. Figure ^ shows the radial dependence 
of the wave with A = Lsobo{r). The phase velocity for the slow 
magnetosonic modes is approximately the velocity of sound, 
so that inward running waves are advected outward in the su¬ 
personic part of the wind. The amplification of the waves is 
strongest close to the star, where the Sobolev length is short. 
The phase velocity of the fast modes goes with for large 
radii since B = ^ r~^. This might lead to a small velocity 

for outward running shocks at large radii. 

Figures show the same plots for model C - with mag¬ 
netic field and rotation. The crucial point is that the magnetic 
field and the amplifying stellar radiation are not parallel any¬ 
more. For large radii they are even perpendicular. Therefore 
the fast magnetosonic wave modes, which are most interesting 
for us, are amplified as well. Figure |[b shows six modes with 
different phase velocities. Two of them, the Alfvenic modes, 
show no dependence on wavelength. They are unaffected by 
the radiation field and therefore stable. Figure ||.a shows the 
amplification timescales for the magnetosonic waves. The fast 
magnetosonic waves grow approximately one order of magni¬ 
tude faster than the slow waves. We can therefore expect that 
these waves will dominate. The crucial point is that they are 
much faster than the unperturbed wind - especially close to the 
star. Inward running fast waves will therefore not be advected 
away from the star. Figure ^b shows that the phase velocity for 
the fast magnetosonic waves remain high for large radii. This 
velocity is a lower lim it for the velocity of outward running 
shocks. Rybicki et al. ( 1990 ) showed that nonradial perturba¬ 
tions in the stellar wind are damped close to the wind. But the 
magnetic field of our model C is mostly tangential already at 
the stellar surface. For the fast magnetosonic modes \Sv^/Svr\ 
is 0.5 at the stellar surface and 0.35 at r = 2R. We expect there¬ 
fore, that the effects found by Rybicki et al. will influence but 
not completely dampen the magnetosonic waves. We emphasize 
that fast magnetosonic waves propagate fastest perpendicular to 
the magnetic field with uph = {v\ -f 

It is very speculative to draw conclusions for nonlinear 
waves and shocks from a linear stability analysis. But we 
showed that a magnetic field has a strong influence on the 
linear analysis. Therefore we expect a strong influence of the 
magnetic field on nonlinear perturbations as well. Our specula¬ 
tive scenario for waves in the wind of hot stars is the following: 
Waves are predominantly generated close to the star where the 
stellar radiation field is strong and the Sobolev length is short. 
Waves with a wavelength short compared to the Sobolev length 
are generated predominantly, because they have the shortest 
amplification timescale. Inward running waves run into the star 
and disappear because their phase velocity is higher than the un¬ 
perturbed wind velocity. Outward running waves, which have, 
for short wavelengths, the same growth timescale than inward 
running waves, can run over many stellar radii and grow. For the 
outward running waves Sp and are in phase. Therefore they 
can steepen into forward shocks with regions of high density at 
high velocity and influence the measurements for the terminal 
velocity and mass loss rate. 
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log(/l/ Lsobil Lsoijt) 

Fig. 1. Amplification timescales (a) and phase velocities (in the frame of the unperturbed wind) (b) versus wavelength for model A&B and 
r = 2R. In model A only modes 3&4 exist. Modes missing in Fig. a) are stable. Mode 3 with A ~ Lsobo has the shortest amplification timescale 
and therefore will dominate the wind. These inward running waves are advected outward with Vr = 811km s“* and steepen into reverse shocks. 



0 10 20 30 40 50 0 10 20 30 40 50 


r[Rl i-LRj 

Fig. 2. Amplification timescales (a) and phase velocities (in the frame of the unperturbed wind) (b) versus radius for model A&B and A = Tsobo. 
In model A only modes 3&4 exist. Modes missing in Fig. a) are stable. The wave amplification is strongest close to the star, where the radiation 
field and wind acceleration are strong. Inward running waves originating there will be advected outward and steepen into reverse shocks. 
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Fig. 3. Amplification timescales (a) and phase velocities (in the frame of the unperturbed wind) (b) versus wavelength for model C and r = 2R. 
Modes missing in Fig. a) are stable. In the long wavelength limit all modes are stable. In the short wavelength wavelength limit the fast 
magnetosonic modes (1&6) grow approximately one order of magnitude faster as the slow magnetosonic modes (3&4). We expect that fast 
magnetosonic waves in both directions are equally dominant in the wind of model C. Since the magnetosonic waves are much faster than the 
unperturbed wind (Vr = 294km s“'), the inward running waves will not be advected outward. The stable Alfvenic modes (2&5) show no 
dependence on wavelength. 




Fig. 4. Amplification timescales (a) and phase velocities (in the frame of the unperturbed wind) (b) versus radius for model C and A = O.lLsobo. 
The Alfvenic modes (2&5) missing in Fig. a) are stable. The wave amplification is strongest close to the star, where radiation field and wind 
acceleration are strong. Fast magnetosonic waves originating here will dominate the wind. They have the shortest amplification timescale for all 
radii and a high phase velocity even at large radii. This may lead to fast shocks running outward, which have a strong influence on observation. 
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4.2. Nonradial Waves 

The picture drawn in the last section does not change signifi¬ 
cantly when nonradial waves are taken into account. There are 
two new effects compared to the radial case. A third wave mode 
appears in model A as predicted by Eq. ^ This inward running 
mode has a long amplification timescale and a relatively low 
phase velocity. Therefore we do not expect a significant contri¬ 
bution of this mode to the situation in the wind of model A. For 
model B we find that Alfven and fast magnetosonic waves have 
different phase velocities now. The fast magnetosonic waves are 
now amplified, too. But their amplification timescale is much 
longer than for slow magnetosonic waves. 

5. Consequences from our model 



Fig. 5. Expansion time versus radius. Model C has a slower accelera¬ 
tion close to the star and a lower terminal velocity. A fast magnetosonic 


wave with an amplification timescale of about 2500s (cf. Fig ^a) grows 
by a factor of « e"^®, while the wind expands from r = li?tor = 1.5R. 
The electron scattering opacity is at r = 1.3R. Therefore we can 
expect shocks alrea dy at t his radius, where line observation starts. The 
objections of Lucy (1984) do not change this situation qualitatively. 


In this paper we do a linear stability analysis and show that the 
waves in the stellar wind can help to understand the observa¬ 
tions. In order to calculate quantitative results, which can be 
compared with observations, it would be necessary to analyze 
the detailed properties of the shocks resulting from the waves 
found in this paper. But we can derive some estimates from our 
calculations. 

From Fig. Qb we see that the phase velocity of the outward 
running waves and possibly shocks in the rest frame of the star 


is at least twice the terminal velocity of the unperturbed wind. 
In the outward running waves the oscillations of Vr and p are 
in phase. Therefore a significant amount of matter, but presum¬ 
ably not all matter will escape at this or a higher velocity, if 
the outward running waves steepen into forward shocks. The 
terminal velocity will then be overestimated at least by a fac¬ 
tor of two in the observation. For our model C this would be 


t'oo.ob.s ~ 1500km s~ ^. Krolik & Raymond ( 1985 ) found that in 
a nonmagnetic wind shock shells are running much faster than 
the phase velocity of the waves. In an unperturbed magnetic 
wind model such a high value for Voo combined with an rea¬ 
sonable high value for M can only be obtained with a very high 
magnetic field and fast rotation, which leads to a spin down 
problem. 

The influence of our model on M is more difficult to esti¬ 
mate. We gain a real factor on 28 in M between our model A&B 
and model C even in the unperturbed wind due to the driving 
force of the rotating magnetic field. Furthermore the observa¬ 
tion of M is influenced by the clumping of the wind matter. But 
to calculate the clumping factor <p^> / < p at large radii, 
where radio observations are made, it would be necessary to do 
a nonlinear calculation including large radii, because the waves 
steepen into shocks very rapidly due to the short amplification 
timescales. This is beyond the linear model presented here. 

Magnetic rotator models for hot stars are often criticized 
for their short spin down times. We reduced this problem by 
using the equations of Biermann & Cassinelli ( 1993| ) for our 
unperturbed wind model C. These equations do not assume that 
the magnetic field is radial close to the star. Therefore we find 
a smaller value for the Alfvenic radius Ra, which controls as 
a lever arm the angular momentum loss. The magnetic field 
increases the angular momentum loss in the equatorial plane 
compared to a nonmagnetic star with the same rotation rate by 
(Ra/R)^, which is 3.2 for our model C. 

In order to decide, wheter the increased angular momentum 
loss would rule out a rotating magnetic field as an important 
component of the wind of these stars, one would need a detailed 
stellar evolution model. Such a model must e.g. consider the an¬ 
gular momentum structure and evolution inside the star, which 
is yet unknown in the case of a strong magnetic field. It is also 
not clear, wheter the magnetic field is present in the wind from 
the beginning of the stellar evolution, or wheter it is generated 


later, as was suggested by Biermann & Cassinelli (19931 and 


Biermann (1997) even for stars with convective cores and ra¬ 


diative envelopes. The angular momentum loss from the wind 
outside the equatorial plane is also not calculated yet. Due to 
these serious uncertainties in this issue most people (including 
the authors) us e the very simple spin down formula of Friend 
& MacGregor (|1984 Eq. 38) 


J 3 R^ M 3 i?2 
TJ --- 7 — 


SR^M 


(24) 


to roughly guess the spin down timescale of the star. This for¬ 
mula guesses the angular momentum of the star by assuming, 
that the star is a rigidly rotating sphere of constant density. All 
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these three assumptions are of course wrong. The angular mo¬ 
mentum loss from outside the equatorial plane is just guessed 
from the value in the equatorial plane using a factor of 2/3. 
Even if we could accurately calculate the current spin down 
J/ j and mass loss M/M timescales, tj < tm does not au¬ 
tomatically mean that post main sequence star does not rotate 
significantly. E.g. an artificial nonmagnetic star {R — Ra), 
wich fullfills the assumptions of Eq. mentioned above , has 
no internal angular momentum transport and a constant density 
over this lifetime, has tj — 3/5tm < tm, but the rotation 
rate fi remains constant over the whole lifetime of such a star. 
The reason is that the outer layers of the star have, due their 
larger average distance from the rotation axis, more angular 
momentum per unit mass than inner layers in spite of the rigid 
rotation. This example shows that Eq. is just not accurate 
enough to verify or rule out a model with T, 7 (Eq. « tm 

within an order of magnitude. Eor our model C we find using 
Eq. = 2.5 X lO^yr and tm = 1-4 x lO^yr. 


6. Conclusions 


In this paper we analyzed the interaction between a magnetic 
field and linear waves induced by the radiative instability. We 
found both models complement each other. The magnetic field 
suppresses the inward running waves, which dominate in non¬ 
magnetic winds. This may allow the outward running w aves to 
support the unperturbed wind as described by Koninx ( 1992 1 
and to form high density shock shells running out at a high 
speed. These shock shells may explain the high terminal veloc¬ 
ities measured in winds of massive stars. The outward running 
waves will also lead to an overestimation of M due to wind 
clumping and the overestimated Voo- Wind clumping also oc¬ 
curs without a magnetic field. But in this case the resulting 
shock shells will run inward; and the argument about the over¬ 
estimated Voo would not apply. The overestimated M and Voo 
put unnecessarily strong restrictions on fast magnetic rotator 
wind models. We argued that even for a magnetic field with 
Bro = 500G the spin down time is consistent with the life¬ 
time of the star inferred from the mass loss rate considering the 
uncertainties in the stellar structure. Erom the observation of 
nonthermal radio emission in many OB and Wolf-Rayet stars 
we know that these stars have a nonneglible magnetic field. Eur- 
ther direct observations are necessary to infer the actual strength 
of these fields. 

We showed that a luminous fast magnetic rotator model plus 
wind perturbations by waves or shocks can help to explain the 
observed high values for M and Voo without being ruled out 
by the spin down problem. Eurther observation of the magnetic 
field and further theoretical work on the evolution of stellar 
rotation are necessary to evaluate to role of magnetic fields in 
winds of massive hot stars. 
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